Quantum black holes with charge, colour, and spin 
at the LHC 



On 
O 



- 1—1 

X 



Douglas M. Gingricb[j] 



Centre for Particle Physics, Department of Physics, University of Alberta, 
<N ; Edmonton, AB T6G 2G7 Canada 

^ , E-mail: gingrich@ualberta.ca 

m 



Abstract. In low-scale gravity scenarios, quantum black holes could be produced 
at the Large Hadron Collider (LHC) provided the Planck scale is not higher than 
a few TeV. Based on fundamental principles and a few basic assumptions, we have 



O i 1 constructed a model for quantum black hole production and decay in proton-proton 



collisions. By performing a particle-level simulation at LHC energies, we have 
' estimated cross sections and branching fractions for final-state particle topologies that 

i-C • would uniquely identify quantum black holes in LHC detectors. If the Planck scale is 

about a TeV, even with the most pessimistic assumptions, the rates for quantum black 
hole production are estimated to be above backgrounds, and some of the final-particle 
states are not found in Standard Model processes. Our results could form the starting 
| point for a detailed investigation of quantum black holes by the LHC experiments. 
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1. Introduction 

Models of low-scale gravity [H [21 [31 HJ |5j [6] allow for the production of small black holes 
in particle collisions El E]- The available energy must be well above the fundamental 
Planck scale for the semiclassical description of black hole production and decay to 
be valid. Based on current experimental and phenomenological limits on the Planck 
scale pm HH H21 H31 OS HS1 HSl EZl HE! HSl 12m I2H 1221 1231 I2^l I2S], it is unlikely that 
semiclassical black holes will be accessible at energies produced by the Large Hadron 
Collider (LHC). However, if the Planck scale is low enough, quantum black holes may 
be produced in abundance at the LHC [261 [271 [28] . 

In spite of not having a complete theory of quantum gravity, it is still possible 
to gain insight into the signatures of quantum black holes at the LHC based on some 
fundamental principles and a few assumptions [29]. The current thinking is that a black 
hole is formed if two partons from a proton-proton collision at the LHC satisfy the hoop 
conjecture [30]. It is natural for the local gauge symmetries of QCD colour and electric 
charge to be conserved during the black hole formation and decay processes [28]. The 
black hole thus inherits the colour, charge, and angular momentum from the parton pair 
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that formed it. Hence, black holes can be characterized by QCD colour, electric charge, 
and angular momentum at the LHC. 

Quantum black holes are expected to differ from semiclassical black holes in a 
number of ways. Semiclassical black holes are thought to decay thermally at the Hawking 
temperature and be well described by black hole thermodynamics. In this canonical 
formalism, the back-reaction on the spacetime during decay is negligible. However, 
as the black hole mass approaches the Planck scale, the back-reaction will become 
significant and the black hole can no longer be considered in thermal equilibrium with 
its radiation. Microcanonical corrections can help improve the decay description until 
quantum mechanical effects take over. When the Compton wavelength of the black 
hole becomes larger than the Schwarzschild radius, quantum effects may impart some 
particle characteristics onto the black hole. At this point in the decay it is difficult 
to view the black hole as still having a well defined temperature and any significant 
entropy. Thus it is unlikely that quantum black holes produced near threshold will 
decay thermally [3U E21 EE] . 

We will consider black holes produced with mass just above the Planck scale. In 
this regime, it is unlikely that the decays will be thermal, and we might expect quantum 
black hole decays into only a few-particle final states to dominate. The formation and 
decay will occur over a small region of spacetime, and the quantum black hole might 
be viewed as a strongly coupled resonance or a gravitationally bound state. Since the 
production and decay are short- distance processes, and the quantum black hole is likely 
to have colour, the QCD hadronization process would occur after the decay. 

With this simplified picture of quantum black hole production and decay, we 
perform a study of the decay signatures that might be expected at the LHC. Our work 
builds on Ref. [291 EE] where these ideas, to our knowledge, were first discussed. We 
explicitly consider only ADD type black holes. However, there is a range of mass scales 
for which almost flat five-dimensional space is an applicable metric for Randall-Sundrum 
type-1 black holes. Some novel signatures that do not often occur in other beyond the 
Standard Model physics scenarios will be examined. Another unique feature is that, if 
black holes are produced, these signatures will occur at huge rates when compared to, 
even, Standard Model processes. 

We reminded the reader that there are inevitable limitations in any model of 
quantum black holes since this is precisely the regime in which gravity becomes strongly 
coupled and the theory is no longer perturbative. Thus our seemly predictive results 
should be viewed as a dimensional analysis in which it is hope that some of the 
extrapolations from the classical domain will carry over to quantum black holes. With 
the startup of the LHC, such phenomenological studies are of interest and of value. 

This paper is structured as follows. In Sec. El we consider the electric charge, QCD 
colour, and spin of black holes produced in proton-proton collisions. A model for the 
production of quantum black holes at the LHC is discussed in Sec. |3] and their decays 
in Sec. HI In Sec. El we estimate cross sections and discuss the topologies that might be 
observed in experiments. 
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2. Quantum black hole states 

Quantum black holes can be classified according to their SU(3) C and U(l) em 
representations. Since we are considering proton-proton collisions, the allowed particles 
forming the black hole are quarks, antiquarks, and gluons. Nine possible electric charge 
states can be formed: ±4/3, ±1, ±2/3, ±1/3, 0. The 4/3 charge state can only be formed 
by quark pairs. The 2/3 charge state can be formed either by an antiquark-antiquark 
pair or a quark-gluon combination. The 1/3 charge state can be formed either by a 
quark-quark pair or an antiquark-gluon combination. The 1 charge state can only be 
formed by a quark- ant iquark pair, and the charge state can be formed by a quark- 
antiquark pair or a gluon-gluon pair. The composition of the black hole negative charge 
states can be enumerated similarly. 

The possible colour states of two partons are 



3©3=8©1 (1) 

3©3=6©3 (2) 

3©3=6©3 (3) 

3 © 8 = 3 © 6 © 15 (4) 

3©8 = 3©6ffil5 (5) 

8©8 = 1 S ©8 S ©8 A © 10©T0 A ©27 5 . (6) 



Since black holes form representations of SU(3) C , they are predominantly coloured, 
but can occur as colourless singlets. Table [H lists the possible charge and colour state 
combinations that could be produced. These states are not unique to quantum black 
holes but also apply to semiclassical black holes formed by two-parton collisions. 

Table 1. Possible quantum black hole states of different electric charge (superscript) 
and colour representation (subscript). For the electric charge states, there are a 
corresponding 10 states of opposite electric charge and "opposite" colour representation 
which are not shown. 

qbhI /3 qbh* /3 

QBH{ QBHl 

QBH^ 3 QBH^' 3 QBH^ 3 

QBH^ 3 QBHl /3 QBH^ 3 

QBH° QBHl QBH% Q BH jo Q BH 27 

At energies of the fundamental Planck scale Mn, the size in spacetime of the 
incoming partons and the gravitational radius r g of the black hole are both of order 
M^ 1 . Since the impact parameter is of the same size as the gravitational radius, the 
angular momentum of the two-particle system, J < Mr g , would be of order unity. Thus, 
a possible semiclassical black hole spin-down process is unlikely to apply to quantum 
black holes. This statement relies on the hoop conjecture of classical gravity being 
applicable in the strong gravity regime. 
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For simplicity, we will consider the initial angular momentum of the quantum black 
hole to be due entirely to the spin states of the incoming partons, and ignore the 
possibility of an initial small orbital angular momentum due to an impact parameter. 
The procedure outline here can be used to include black holes with orbital angular 
momentum in the future, if desired. When considering spin states, we will assume 
massless partons. Since we are ignoring angular momentum (transverse to the helicity 
axis), the quantum black hole spin will be parallel to the parton helicity axis, or the 
beam axis at the LHC. This is quite different from semiclassical black holes with angular 
momentum. In this case, the convention is to add the small amount of spin angular 
momentum to the normally larger orbital angular momentum, which is transverse to 
the helicity axis [33]. The quark-quark states can form spin-0 and spin-1 quantum 
black holes. The spin-1 state is three times more likely to form than the spin-0 state. 
The quark-gluon states can form spin-3/2 and spin- 1/2 quantum black holes. However, 
because the partons are massless, not all spin combinations are possible. The spin- 
3/2 state is twice as likely to form as the spin- 1/2 state. The gluon-gluon states can 
form spin-0, spin-1, and spin-2 quantum black holes. Again because the partons are 
massless, not all spin combinations are possible. There are no ±1 states along the 
spin axis in the spin-2 case. The relative spin-2, spin-1, and spin-0 probabilities in 
gluon-gluon collisions are 7:3:2. It may well be that charge, colour, and spin are not 
totally independent. The over all wave function of the two parton system may have 
to be symmetric or antisymmetric under interchange of two identical quarks or gluon, 
respectively. 

3. Production of quantum black holes 

The cross section for black hole production is not known. Based on classical arguments 
and only one available scale, the cross section is most often taken to be the geometrical 
cross section a ~ nr^, where r g is the gravitational radius of the two-particle system (see 
Ref. [31] for a review). This cross section assumes a radiationless production process in 
which all the energy of the two-particle system goes into forming the black hole. Various 
calculations have been performed to predict the amount of gravitational radiation in 
the production process, but for higher- dimensional gravity only the trapped surface 
approach has yielded numerical results [351 EHJ E7]. The effect of radiation on the cross 
section is smallest for low black hole masses [31]. For the quantum black holes that we 
will consider, the trapped surface cross section is about 50 time less than the geometrical 
cross section. However, for higher- dimensions and non-zero impact parameter, these 
calculations can only be considered as lower bounds on the cross section. In addition, 
the calculations have ignored parton charge, parton spin, and parton finite size (see 
Ref. [381 EH1 HQl HI] for progress in these areas) . The calculations are based on classical 
arguments, and since it is not known how applicable these might be in the quantum 
regime, we will not consider them except for to realized that the cross sections might 
be lower by as much as about two orders of magnitude. 
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Arguments have been made for why the geometrical cross section, to within a 
small numerical coefficient, should be applicable to quantum black holes [121 H3] . 
Alternatively, if one embeds a quantum theory of gravity into string theory, the cross 
section is reduced due to the quantum size of the string [JT]. The energy threshold 
for quantum black hole production would be somewhat larger than Mo by a factor 
that depends on the ratio of string length to Planck length, or the string coupling in 
weakly-coupled string theory. The amount of reduction in the cross section will remain 
unknown until these fundamental constants have been determined. 

We assume the production cross section for quantum black holes can be 
extrapolated from the cross section for semiclassical black holes. It is a challenge for 
the experiments to measure or set a limit on the black hole cross section. This might 
be one of the few numerical quantities the experiments can address in an almost model 
independent way. 

The gravitational radius r g of a quantum black hole of mass M becomes 



where D is the total number of spacetime dimensions, and k(D) is a numerical coefficient 
depending only on the number of dimensions and the definition of the fundamental 
Planck scale; the PDG definition of the Planck scale is used in this study: 



We now address the question of over what mass range is a black hole a quantum 
object? It is common to take the validity of the semiclassical black hole to be when the 
Compton wavelength of the colliding particles lies within the gravitational radius. We 
turn this condition around and consider it to be the upper mass bound for considering 
the black hole to be in the quantum regime. With this mass restriction, we stay away 
from the semiclassical regime, where black hole thermodynamics and thermal decays 
occur. This bound also ensure that the initial angular momentum of the quantum black 
hole is close to unity. If we do not take the upper mass requirement into account, the 
tails of distributions (particle momentum, for example) might be artificially altered due 
to semiclassical black hole decays. We consider the lower mass, or threshold, to be 
when the mass of the black hole is equal to the inverse of its radius. Using the PDG 
convention for the Planck scale gives 



We notice that the minimum mass is below the fundamental Planck scale. This is an 
artifact of the definition of the Planck scale. Using the more intuitive Dimopoulos- 
Landsberg definition of the Planck scale gives a minimum mass always above the 
fundamental Planck scale. On the other hand, the accelerator experiments have set 




(7) 




(8) 




(9) 
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limits on M D of M D > 1 TeV, and hence a quantum back hole would be require to have 
mass above only about 500 GeV. If Mjj is about 1 TeV, could it be that the Tevatron 
has failed to observe the effects of relatively low-mass black holes? Assuming this is 
not the case, we will take the threshold for quantum black hole production to be the 
Planck scale in the PDG definition, and consider a fixed mass range for all dimensions: 
M min = Mfl and M max = 3Mp . The maximum mass is lower than that indicated by 
Eq. (Q, but it is the lowest common choice where semiclassical approximations are 
expected to break down, and thus represents a very conservative maximum mass for our 
model. 

Only a fraction of the total centre-of-mass energy y/s in a proton-proton collision 
is available in the hard scattering process. We define s, where x a and 

Xb are the fractional energies of the two partons relative to the proton energies. The full 
particle-level cross section a is given by 




where a and b are the parton types in the two protons, and f a and fb are the parton 
distribution functions (PDFs) for the proton. The sum is over all the possible quark 
and gluon pairings that can make a particular quantum black hole state. The parton 
distributions fall rapidly at high relative energies, and so the particle-level cross section 
also falls at high energies. 

When describing quantum black hole production, it will be more informative to 
specify the two partons pi and p 2 that went into the formation of the black hole, and drop 
the colour representation specifier by summing over the possible colour representations. 
This makes sense since the partons hadronize after the quantum black hole decays. The 
relative probabilities of each colour representation will be accounted for in the decay 
branching fractions. The resulting 14 different processes can be used to obtained any 
desired inclusive result by multiplying the cross section by the branching fraction and 
summing the relevant processes. 

To obtain numerical results, we have used the parameters shown in Table EJ You 
may assume these parameters were used in calculations unless told otherwise. The 
CTEQ6L1 parton distribution functions jB] were used with a QCD scale of Q = l/r g . 
The cross section spectrum of quantum black hole states is shown in Fig. [TJ It is 
interesting to note that the highest cross section is given by u-g collisions. This shows 
the importance of the gluon contribution at the parton kinematics and QCD scale that 
we are using. As expected, the inclusive cross section, of about 130 nb, is significantly 
higher than the semiclassical black hole cross section. This value should be considered 
as an upper limit. 

The cross section is sensitive to the choice of PDFs and QCD scale. Table [3] shows 
the inclusive quantum black hole cross section for four different choices of PDFs and 
two different QCD scales. A difference of about 7% is observed for the different PDFs. 
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Figure 1. Cross section spectrum for quantum black holes of different electric charge 
and parton progenitors, u stands for all up- type quarks (u, c,t), d stands for all down- 
type quarks (d, s, b), and similarly for the antiquarks. 



Our calculations use Q = l/r g for the QCD scale. If Q = M is used, the cross sections 
are about 13% lower. 

Figures [2] shows how the total inclusive cross section changes with different number 
of dimensions and different values for the Planck scale. If the LHC detectors can 
collect a few hundred pb _1 of data in the first years, they will be able to produce 
significant numbers of quantum black holes for Planck scales as high as about 5 TeV, 
if the geometrical cross section is valid. Also shown in Fig. |2] are the trapped surface 
lower-bound cross sections. These cross sections are about 10 to 10 4 lower than the 
geometrical cross section over the range 1 < Mr, < 4 TeV. In this case, we see that 
there is very little dependence on the number of dimensions. In calculating the trapped 



Table 2. 


Default parameters used in calculations. 


Value 


Symbol 


Description 


10 


D 


total number of dimensions 


1 TeV 


M D 


fundamental Planck scale 






(PDG definition) 


14 TeV 




LHC centre-of-mass energy 


M D 


M min 


minimum quantum black hole mass 


3M D 




maximum quantum black hole mass 
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Table 3. Inclusive quantum black hole cross section in nanobarns for different parton 
density functions and different QCD scales. 



QCD 


Parton Density Functions 




Scale 


CTEQ6L1 CTEQ5L CTEQ5M 


MRST98 


Q = l/r g 


132 136 159 


150 


Q = M 


115 119 139 


132 



surface cross section, the quantum black hole mass has been limited to the range of 
My) < M < 3Me>- That is, if the mass before gravitonal radiation is above 3Mo it is 
not included in the total cross section calculation. 



3 — i — | — i — i — i — i — | — i — i — i — i — I — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — c 




Planck Scale [TeV] 



Figure 2. Total proton-proton cross section for different number of dimensions D and 
different fundamental Planck scales. The solid lines are totally inelastic cross sections 
and the dashed lines are trapped surface cross sections. M max = min(3MD, y/s) has 
been used. 

Because of the falling product of parton distribution functions with parton-parton 
centre-of-mass energy the most probable value for the quantum black hole mass 
is Mb- At this value, the gravitational radius, and hence cross section, is independent 
of the quantum black hole mass. Thus the shape of the cross section at the lowest 
possible masses is almost independent of the parton cross section and is determined 
predominantly by the parton distribution functions. Since each quantum black hole 
state is made from unique valance quarks, sea quarks, or gluons, the product of parton 
density functions can be very different. Figure [3] shows the mass distribution for two 
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rather different quantum black hole states. Although gluons contribute the most at 
low masses, the u-quarks (valance quarks) contribute the most above masses of about 
I.BMd- Over the mass region 1 < M < 3 TeV, the uu — > QBH 4 ^ differential cross 
section falls as a power law in the mass with exponent —2.1, while the gg — > QBH° 
cross section falls as a power law in the mass with exponent —4.6. 




Quantum Black Hole Mass [TeV] 

Figure 3. Quantum black hole mass distribution for two different quantum black hole 
states. 



Since quantum black holes are predominantly produced near threshold, their 
kinematics are different from the kinematics of semiclassical black holes. Figure H] shows 
the kinematics of quantum black holes. The mean mass is about 1.5Md, the mean 
energy is about 2M£>, and the mean momentum is about 0.9M/). The distribution of 
values of the Lorentz variables and 7 can be easily explained. They depend only on 
the x a and Xb values of the two partons forming the black hole, or equivalently on the 
s and s values. (3 and 7 are given by 

(5 = and 7 = . (11) 

•Ea ~T~ %b v ^ a *^b 

For symmetric collisions (x a ~ Xb), (3 ~ and 7 w 1. For highly asymmetric collisions 
(xi w 1 and x 2 ~ x m - m ), 



/5 ^ ^ — !1 ^ 0.99, 7 ^i±i^^L^:%^7. (12) 
s + S 1 2yfsl 2v^ 2 TeV K ' 

Thus, based on Fig.HJwe can see that the collisions are more asymmetric than symmetric. 
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Figure 4. Distributions of some kinematic variables of quantum black holes. 



4. Quantum black hole decay 



Up until now, we have visualized the hard scattering process as occurring through a 
quantum black hole resonance; in analogy to a particle resonance. It is informative to 
estimate the width of this possible resonance. For a well defined resonance, we will use 
the criteria that r 3> 1/M, where M is the initial quantum black hole mass and r is 
its lifetime. It is assumed that the quantum black hole decays entirely without leaving 
a remnant. When estimating the lifetime, it is usual to assume Hawking emission 
of massless particles during the entire decay process. One can start with the power 
spectrum for black hole emission or go directly to the generalized Stefan-Boltzmann 
equation and relate the Hawking temperature to the Schwarzschild radius. The total 
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power emitted will be equal to the rate of decrease of the black hole mass. The resulting 
differential equation is then integrated over the entire lifetime of the black hole and the 
result is given by 



where the form of C depends on how the emissivities of the different fields are normalized. 
In any case, C involves a sum over the particle degrees of freedom for each spin weighted 
by the integrated power emitted into each spin field. The power factors are greybody- 
modified thermal power spectra proposed by Hawking [15] . 

Since we are considering non-thermal quantum black hole decays, the applicability 
of Eq. ffl3|) is questionable. However, we will consider it as an indicative estimate of 
the quantum black hole lifetime. Since quantum black holes are produced with mass 
close to the Planck scale, we expect tM > 1/C. Thus our requirement for a resonance 
translates into C < 1 or large mass (at which point the black hole is semiclassical) . 
The coefficient C increases with the number of particle degrees of freedom, the value of 
the greybody factors, and the number of dimensions. If greybody factors are ignored, 
C ~ 2 for D = 10 and does not drop below unity until D < 8. If we use the greybody 
factors calculated according to Ref. [16] and only consider emission on the brane, C ~ 5 
for D = 10, and C is still greater than two for seven dimensions. If we include the 
emission of gravitons into the bulk using the results of Ref. [IT] , the number of degrees 
of freedom increases, particularly for higher dimensions, and C becomes even larger. 
Thus, unless we push the limits of the quantum black hole mass into the semiclassical 
regime or restrict our considerations to low number of extra dimensions (those likely 
to be excluded by experiments), a quantum black hole probably can not be viewed 
as a particle resonance. Some authors have reached similar conclusions [1H1 EH] . The 
concept of a quantum black hole state still has the useful purpose of labeling the possible 
interactions. 

In the decays of quantum black holes, one expects the number of particles in the 
final state to be small. Using arguments similar to those used to estimate the lifetime, 
the average number of particles emitted from the black hole during Hawking evaporation 
can be estimated to be 



where S is the initial entropy. The form of p depends on how the emissivities of the 
different fields are normalized. If greybody factors are ignored, the average multiplicity is 
between 1.0 to 1.3 for 6 to 10 dimensions. If greybody factors are included and emission 
is restricted to the brane, there is no significant change in the average multiplicities. 
Allowing graviton emission into the bulk gives the numbers shown in the second column 
of Table HI The mean multiplicity increases with the number of dimensions to a 
maximum for nine dimensions and then decreases. This is due to the interplay between 




2(13-2) 



(13) 




(14) 
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the entropy definition (due to the Planck scale definition) and the ratio of the sum of 
fluxes to sum of powers in the coefficient p in Eq. (JHJ). 

Table 4. Mean number of emissions (N) and probability of each number of particles 
in the final state versus the number of dimensions D. 



D (N) Number of Particles 







2 


3 


4 


5 


6 


5 


0.6 


0.74 


0.21 


0.04 


0.01 


0.00 


6 


0.9 


0.61 


0.28 


0.09 


0.02 


0.00 


7 


1.1 


0.54 


0.30 


0.11 


0.03 


0.01 


8 


1.2 


0.51 


0.31 


0.13 


0.04 


0.01 


9 


1.3 


0.50 


0.31 


0.13 


0.04 


0.01 


10 


1.2 


0.51 


0.31 


0.13 


0.04 


0.01 


11 


1.1 


0.54 


0.30 


0.11 


0.03 


0.01 



Fluctuations about the mean multiplicity can be described using a Poisson 
distribution [IHJ [25]. The Poisson distribution can also be used to estimate the 
relative probabilities of two-particle, three-particle, four-particle, etc. final states. When 
calculating the relative probabilities, we have removed the case of zero particle emission 
and renormalized the Poisson distributions. We consider single-particle emission to 
represent two-particle decay in that the remaining black hole is considered to be the 
second particle. Using these concepts, the probabilities for different number of particles 
in the final state are shown in Table HI For D — 7 to 11 there is little dependence on 
the number of dimensions. Approximately 50% of the decays are two-particle, while 
three-particle and four-particle decays are not insignificant. The multiplicities depend 
on the definition of the Planck scale. For the Dimopoulos-Landsberg definition and the 
case of D — 10, the mean multiplicity is 0.4 and the probability of a two-particle decay 
is about 80%. The two-particle state is enhanced because, by using a different definition 
of the Planck scale, we are effectively considering a lower threshold for quantum black 
hole production. In what follows, we will only consider two-particle final states. Clearly, 
an extended analysis should include higher multiplicity decays. 

For the decay kinematics, we have used two-particle phase space decays. The decay 
particle kinematics are shown in Fig. [5] In the Lorentz and 7 distributions, we have 
excluded zero-mass particles and low-mass particles: electrons, muons, u-quarks, and 
d-quarks. Since we are using two-particle decays, the distributions are representative 
of all particle types, except for the distributions of the Lorentz variables /3 and 7. The 
most probable particle energy, momentum, or transverse momentum is about Md/2. 
Most particles will remain in a detector and are highly relativistic. 

In quantum black hole decays, we consider the colour and charge to be conserved, 
but make no similar assumptions about global charges like baryon or lepton number. We 
consider gravity to be totally democratic and couple to all flavours of quarks, leptons, 
and gauge bosons equally. (See Ref. [501 EH E2] for a discussion of non-democracy in 
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Figure 5. Distributions of some kinematic variables for decay particles from quantum 
black holes. Massless particles, electrons, muons, w-quarks, and d-quarks have not 
been included in the Lorentz /3 and 7 distributions. 



a model with many particle species.) Although perhaps debatable, we assume Lorentz 
invariance holds, and only allow decays that conserve total angular momentum. 

Experimental bounds on the effects of higher-dimension operators might limit the 
parameter space over which we can consider decays violating global charges. The limits 
on flavour-changing lepton decays and proton decay should be taken into consideration. 
Unfortunately, it is not clear what the relationship is between the Planck scale and 
the scale occurring in higher-dimension operators. Arguments have been made for why 
the experimental constraints are serious but do not necessarily rule out the allowed 
parameter space that we will consider [29j [28]. It is also possible that the global 
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symmetries may be gauged. 

Baryon and lepton violating processes may also be mediated by virtual black 
holes and thus affect higher order diagrams [53]. The probability of proton decay 
mediated by a virtual black hole is very small for regular high-scale gravity in four 
dimensions [54J. However, in low-scale gravity, the baryon number violating processes 
through formation of an intermediate virtual black hole can become sizable. Various 
methods have been proposed to suppress these baryon and lepton violating processes to 
agree with experimental limits [551 156] . 

Quantum black holes have colour, and size of about M D X < 1 TeV" 1 . The decay 
products can be partons in coloured representations, which can travel a distance of 
about the QCD scale of a Fermi before they have to hadronize to form colour singlets. 
Thus, we only allow parton decay final states that can have the colour representation of 
the quantum black hole, or equivalently, the colour representation of the initial parton 
progenitors. Since it would be difficult to identify the quantum black hole colour in 
experiments, we will consider the different singlet and octet representations (see Eq. (CQ) 
and 06])) to be the same when counting the number of degrees of freedom. 

The particle content of low-scale gravity is uncertain. In fact, a large number of 
dark sector species may reveal themselves in the strong gravity regime J5j (6J [52]. A 
Higgs boson might be discovered and there is likely to be a graviton. In addition, 
how the neutrino sector couples to gravity is an interesting topic. Now that at least 
one neutrino has mass, it is possible that gravity will couple to neutrinos of both left- 
and right-handedness equally. In addition, whether neutrinos are Dirac or Majorana 
particles will change the number of available final states. However, the biggest factor 
effecting the final state particles is if the global symmetries of the Standard Model remain 
good symmetries at the strong gravity scale. In presenting our results, we most often 
consider two extreme particle content models. Our base model allows global symmetries 
to be violated, includes a Higgs and graviton, and assumes the neutrinos are Majorana 
particles of both handedness. Our comparison model is the currently observed Standard 
Model: global symmetries are conserved, there is no Higgs or graviton, and the neutrinos 
are left-handed Dirac particles. When considering the different particle content models 
it is important to realized that coupling per degree of freedom in Hawking radiation does 
not apply in the quantum gravity cases. For quantum gravity, we are simply assuming 
the coupling to all quark and lepton flavours are equal. 

We determined the branching fractions for each quantum black hole state as follows. 
For each state, we wrote down all the two-particle final states that could conserve colour, 
charge, and angular momentum. Some decay states can not conserve all the angular 
momentum modes and were thus weighted accordingly. For each decay mode, we formed 
a product of weights given by the number of members in the colour representation, times 
the number of flavour combinations, times the spin degrees of freedom. The decay 
branching fractions are shown Table [5j 

We reminded the reader that the quantitative results presented in Table [5] and the 
following section are of an illustrative nature. They are based on many assumptions 
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Table 5. Quantum black hole decay branching fractions. BR is the branching fraction, 
C means global symmetries are conserved and V means they may be violated. 
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and extrapolations from the semi-classical regime. 
5. Results 

Using the previously developed model, we have simulated the production and decay 
of quantum black hole events, with each state weighted by its cross section, using a 
new Monte Carlo event generator described in Ref. [57]. To out knowledge no such 
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study has yet been performed. The generators Charybdis [33] and BlackMax [53] 
both include some facility to generate quantum black holes. BlackMax implements 
the suggestions of Meade and Randall verbatim [29]. The decay species are chosen 
according to thermal greybody-modified distributions, which we have argued might not 
be applicable to quantum black hole decays. Charybdis allows the user to set parameters 
so that effectively a quantum black hole decaying to two-particle final states is generated, 
but does not take colour explicitly into account when determining the two-particle final 
state. In both cases, the generators do not allow for the ability to generate a particular 
quantum black hole. They can thus be made to act as inclusive quantum black hole 
generators. In addition, these generators do not allow baryon number violating processes 
that can be subsequently hadronized. The new generator can be used to select the 
specific quantum black hole state to generate, includes all possible decay modes properly 
weighted by conserving SU(3) C and U(l) em , but can violate baryon number. 

The particle PDG identifier codes for the two decay particles in the hard scattering 
are shown in Fig. [6j We notice the expected dominance of quarks and gluons, and the 
charge and baryon asymmetry from having two protons collide, rather than a proton- 
antiproton collision. Figure [6] is for the observed Standard Model particle content and 
can look significantly different depending on the particle content model. If neutrinos 
are chiral, the number of charged leptons and neutrinos are approximately equal. If 
global symmetries can be violated, the number of charged antileptons and antineutrinos 
increases, while the number of antiquarks decreases. Including the graviton is significant, 
while the Higgs boson has little effect. Figure [7] show the frequency of particle 
identification codes when including all these effects. We notice a significant increase 
in leptons and near equality between the different antiquark flavours. 

We have studied the multiplicities of various particle signatures in a detector: jets, 
electrons, muons, photons, and missing energy. For this calculation, we allowed the t, 
W, Z, and H to decay. We called all quarks, gluons, and tau particles jets. On average 
per event there are 2.4 jets, 0.05 electrons, 0.05 muons, 0.05 photons, and 0.2 particles 
(neutrinos and gravitons) that give missing energy. There is very little variation due to 
changing the particle content model. If global symmetries can be violated, the average 
number of electrons, muons, and particles giving missing energy per event are marginally 
higher. This is because of the domination of jets in both cases. 

It is useful to know the average number of jets in events and how the decays of 
the t-quark, W, and Z effect these numbers: 99% of the events have at least one jet, 
87% of the event at least two jets, and 28% of the events at least three jets. If global 
symmetries are violated, the number of events with two or more jets decreases by about 
10%. Parton hadronization, detector effects, and jet-finding algorithms will cause the 
number of jets per event to be higher. 

To study the topologies of the events, we do not allow the final-state particles to 
decay. The percentage occurrence of each topology is shown in Table |6j As expected, 
di-jet events dominate. Also significant are mono-jet, jet plus t-quark, jet plus gauge- 
boson, and jet plus lepton topologies. Noteworthy is the electron plus muon topology 
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Figure 6. Relative occurrence of decay particle types according to their PDG 
identification code. Observed Standard Model particle content with global symmetries 
conserved. 



at the 0.08% level. Such a signal is not usually produced in other beyond the Standard 
Model physics processes. The cross sections are also shown in Table El along with 
predictions of the Standard Model when using PYTHIA 8 [59, 60J. For the Standard 
Model processes, we have restricted the phase-space by requiring the invariant mass 
to be between 1 TeV and 3 TeV, and the transverse momentum in the rest frame of 
the hard scattering process to be between 500 GeV and 1.5 TeV. In most cases, the 
Standard Model cross sections are more than three orders of magnitude lower. So even 
if the quantum black hole cross sections are at the lower bound given by the trapped 
surface calculation, the number of events produced in each decay topology should be 
greater than the Standard Model background. 

The previous analysis used a Planck scale of 1 TeV, which is a hard lower bound 
for D > 6. If one considers how Kaluza-Klein gravitons would affect supernovae cooling 
and neutron stars, the lower bounds on the Planck scale are higher [T8~| [19]. However, 
it should be realized that the bounds provided by astrophysical arguments contain 
significant uncertainties, such as the compactification moduli. For D < 8, the lower 
bound on the fundamental Planck scale is too high to allow quantum black holes to 
be produced by the LHC. For higher dimensions, the bounds on the Planck scale are 
less stringent. For D = 8, M D > 4 TeV and for D > 8, M D > 1.4 TeV. For D > 8, 
the lower bound on the Planck scale is set by the absence of black holes in neutrino 
cosmic ray showers [25]. However, Auger has yet to observe a single neutrino-induced 
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PDG Particle Identification Code 

Figure 7. Relative occurrence of decay particle types according to their PDG 
identification code. Global symmetries may be violated, neutrinos are chiral and 
Majorana, and a Higgs and graviton are allowed. 



shower. Ref. [25] is based on the ratio of vertical and quasi-horizontal neutrino showers. 
The vertical showers are used to normalize the product of the neutrino flux and their 
interaction cross section. The fact that Auger should have seen a few vertical showers 
by now, but have not seen any yet, relaxes the limits on Mn. In addition, the neutrino 
cosmic ray bounds can be evaded in a model of split fermions (61]. Using the results of 
Fig. El we see that for D = 8, a ~ 20 pb and for D > 8, a > 2 x 10 4 pb. Thus, most 
decay signatures in Table |6] would be observable for D > 8, while probably none would 
be observable for D < 8. 

We comment on the D = 8 little more by giving some order of magnitude 

estimates. We assume the ratio of the number of observed events to the experimental 
acceptance is 10 3 . For di-jet events, the over all acceptance might be high, but some 
background can be anticipated and hence on the order of 100 events will need to be 
observed. While in the case of leptons, the backgrounds might be lower, but the 
acceptance will probably be lower since more stringent requirements might be needed to 
reduce lepton fake rates, and ensure the leptons are not accompanied by additional jets. 
Thus about 80 pb -1 of data would be be required to see an anomaly in the Standard 
Model di-jet events, and about 1 fb _1 to see one in a jet plus lepton signal. About 
50 fb -1 of data would be needed for a definitive statement to be made in the electron 
plus muon channel. 
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Table 6. Percentage occurrence and cross section for each decay topology. Columns 
two and three are for quantum black holes that conserve global symmetries and columns 
four and five are for quantum black holes that may violate global symmetries. The 
last column is for Standard Model QCD and clectrowcak processes. 
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6. Summary 

Since many assumptions have been made in developing our toy model, we conclude 
by summarizing them. The model is based on assuming local charges (colour, electric 
charge, and spin) are conserved in the strong gravity regime of black hole production 
and decay near the Planck scale. However, we have not taken any initial orbital angular 
momentum of the black holes into account, but argued that it would be small. We 
have only examined two-body final states and motivated this choice by extrapolating 
semiclassical results to the Planck scale. We have required Lorentz invariance and have 
not considered a black hole remnant - both arbitrary, but working choices. Conservation 
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of global charges like baryon and lepton number has not been required. Predictions based 
on global charge conservation and non-conservation have been presented. 

Experiments at the LHC are currently preparing to search for the effects of low-scale 
gravity below the Planck scale by searching for gravitational scattering and resonances, 
and above the Planck scale by searching for semiclassical black holes. However, there 
has been little guidance in the literature as to the phenomenology of low-scale gravity at 
the Planck scale. This is because of the quantum and non-perturbative nature of gravity 
in this strong-gravity regime. To help guide experimentalists, we have considered black 
hole production and decay near the Planck scale. Based on fundamental principles 
and some assumptions, we have built a model to estimate cross sections and decay 
topologies, albeit aware of its limitations. Although di-jet decays denominate, they are 
unlikely to be the only decays with a significant rate. The jet plus gauge-boson and jet 
plus t-quark topologies account for about 30% of the decays. In models where global 
symmetries need not be conserved, jets plus leptons and mono-jets are also significant 
decays. If the Planck scale is low enough so that the inclusive cross section is of the order 
of 100 nb, then even the decays with low branching fractions will be readily observable. 
Signatures such as opposite sign electron plus muon, with very little else in the detector, 
will be hard to explain by conventional physics. Observing anomalous rates in several 
of the decay channels could help determine the nature of the new physics and lead to 
the discovery of quantum gravity. 
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